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I. INTRODUCTION 



The path integral is a formulation of quantum mechanics equivalent to the standard formu- 
lations, offering a new way of looking at the subject which is more intuitive than the usual 
approaches. Application of the path integrals are as vast as quantum mechanics itself, in- 
cluding the quantum mechanics of a single particle, statistical mechanics, condensed matter 
physics and quantum field theory pj, |2j|. Quantum field theory is the quantum mechanics 
of continuous systems and fully developed in quantum electrodynamics which is the most 
successful theory in physics. The path integral in quantum field theory, both relativistic 
and non-relativistic, play a much more important role, for several reasons. They provide 
a relativity easy road to quantization and to expressions for Green's functions, which are 
closely related with many physical quantities. The close relation between statistical mechan- 
ics and quantum field theory is plainly visible via path integrals. Usually we are interested 
in a quantum field which has to be considered in the presence of a matter field described by 
some bosonic fields. For example, in quantum optics there are some problems where we need 
to quantize the electromagnetic field in the presence of some dielectrics or in calculating the 
effects of matter fields on Casimir forces pHH]. For a consistent quantization of such fields 
we should include the matter fields in the process of quantization. But the main problem is 
that there are few situations where the interested physical quantities, like for example, the 
Casimir force, can be determined in a closed form and finding an effective approximation 
method is quite necessary [7]. 

The Casimir energy is the difference between the energy of the fluctuating field when 
then objects are present and when the objects are removed to infinity. The advent of 
precision experimental measurements of Casimir forces P4T2] and the possibility that they 
can be applied to nanoscale electromechanical devices [131 113] has stimulated interest in 
developing a practical way to calculate the dependence of Casimir energies on the shapes of 
the objects. Many geometries have been analyzed over the years, but the case of compact 
objects has proved rather difficult. One way to calculate the Casimir effect is using the 
multiple scattering formalism, which dates back at least to the 1950. Recently Milton and 
et al., have noticed that the multiple scattering method can yield exact, closed form results 
for bodies that are weakly coupled to the quantum field [I5l IT6] . 

In the present work, we use path integrals to calculate the Casimir force between arbitrary 
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shaped objects foliated in a scalar or vector fluctuating fleld. For this purpose, we start from 
a Lagrangian and develop a perturbative approach to calculating the Casimir force in the 
case of scalar and vector flelds by taking into account the susceptibility function of the 
medium as the expansion parameter. The covariant formulation of the method and also 
its flnite-temperature correction is given which may have applications in dynamical Casimir 
effects fr7\ [18] . Finally, it is shown that the results are consistent with the weak coupling 
limit in multiple scattering formalism as expected [151 UHl [T914M] . The present work is in 
fact a generalization of the previous work |27] as follows: (i) The massless scalar fleld in 
1 + 1-space time is fully investigated and an expression for Casimir force between intervals 
is derived, (ii) The Proca electromagnetic fleld is investigated and an expression for the 
free energy in zero and flnite temperature is obtained for two arbitrarily shaped dielectrics, 
(iii) the approximations are improved and the results are compared with the weak coupling 
regime in the scattering approach. 



II. MASSLESS SCALAR FIELD 



The Lagrangian of a massless Klein-Gordon fleld in iV + 1-dimensional space-time {x = 
(x, t) G ]R^+^) is the simplest fleld theory given by 

= ^dM^)dM^)- (1) 

As usual, let the medium be modeled by a continuum of harmonic oscillators, namely the 
Hopfleld model [25] as follows 



oo 








This modeling can be easily generalized to a covariant form but let us flrst discuss the non 
relativistic case. The interaction between the scaler fleld and its medium is assumed to be 
linear and described by 

oo 

Cint = j dujf{uj,x)Y^{x)Lp{x), (3) 



where /(w,x) is the coupling function between the scaler fleld and its medium. If the 
medium has fllled a flnite space Q, then for (x ^ there is no matter to be coupled to the 
fleld so we should set /(a;,x) = 0. Having the total Lagrangian, we can quantize the total 
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system using path-integral techniques [26] . An important quantity in any field theory is the 
generating functional from which n-point correlation functions can be obtained by taking 
successive functional derivatives. Here our purpose is to find two-point correlation functions 
or Green's function in terms of the susceptibility of the medium. For this purpose the free 
generating-functional can be written as [27] 



2h 



yjd":x.jdtj d"x' / dt'J{x)G'^{x-x')J(x') 



-^Jd'^^Jdifd"^' J dt' f dLoJ^ (xOGO {x-x')J^ {x') 
X e . (4) 

From the free generating-functional we can obtain the interacting generating-functional from 
the well known formula 



W[J,J^] = e Wo[J,J^] 

oo 

X e'^^'^"''^'^^^ ^ dt'J(x)G^(x-x')J{x') 

J d"x fdtj d"x' J dt' J di^Ju,{x)G%(x-x')Juj{x') 

X e . (5) 



where 



f Hkn piko(t-t') 

G»(.-.')-M«(x-x')/^^^^^,— ^, (6) 

Now having the generating functional, the two-point function or Green's function can be 
obtained as 

By using Eq. (j5]) and after some straightforward calculations, we find the following expansion 
for Green's function in frequency variable 

G(x-x',u;) =G°(x-x',w)+ [ ^ziG^x - zi,u)[uj^x{(^,zi)]G\zi - x' ,uj) + 

Jn 

rf"zici"z2G'°(x -zuoo) [oo^xioo, zi)]G'°(zi - Z2, oo) [uj^xi^, Z2)]G°(Z2 - x', u;) + ■ ■ ■ , (8) 




n Jn 



where x(ci;,x) is the susceptibility function of the medium in frequency variable. 
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By using Euler lagrange equations we find equation of motion for the fields and in par- 
ticular we find a Langevin type equation for the scalar field. Specially we can show that the 
Green's function of the system satisfies 



t 



9^9'^G(x-x',t-t') + ^ j rft"x(t-t",x)J^G(x-x',t"-t') = 5(x-x',t-0. (9) 

— oo 

We can also show that the coupling function between the scalar field and its medium is 
proportional to the imaginary part of the susceptibility which is responsible for dissipation 
in the system 

/2(a;,x) = -53[x(u;,x)]. (10) 

71 

The memory function in a homogeneous medium is position independent and therefore, 
Eq.([9]) can be solved easily in reciprocal space 

^(k, u:) = p / . = \,(.) , (11) 

k — w' — u^xi'^) k — u^^^ 

where = 1 + x(w) can be considered as the dielectric function of the medium. From 



Eq.(ll) it is clear that Green's function in the presence of a homogeneous medium can be 
obtained from Green's function of the free space simply by substituting o;^ with o;^^^. Now 



Eq.([9]) in frequency-space can be written as 

[V2_a;2^i^^]G(x-x',u;) = 5(x-x'). (12) 

^0 

In some simple geometries where dielectric function is a piecewise homogeneous function, 
for example for an array of parallel dielectric slabs, Eq.([9]) has an exact solution, see for 
example The function e{uj,x.) in a general piecewise homogeneous medium is defined by 

e(^,x) ^ I e(w), X G 
eo I 1, X ^ 

where fl is the union of regions filled with a matter with the dielectric function e{u). In 



general, finding an exact solution of Eq. ( 12 ) is quite complicated or even impossible and a 
series solution will be useful. 
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III. COVARIANT FORMULATION 



For obtaining a covariant formulation we model the medium by a continuum of Klein-Gordon 
fields and write the total Lagrangian density as follows 

where 



1 

— ( 

2 

oo 

1 



Cres = duid^YMd^yM-u'Y^ix)), 


oo 

dnt = j dujf''{uj,x)Y^{x)d^(^{x), (15) 







are the Lagrangian densities of the scalar field, its medium and their interaction, respectively. 
Here x = (x°,x) belongs to the Minkowski space-time M^^-^. The fields (p{x) and 
are scaler fields and the coupling function f^{uj,x) is a vector field. From Euler-Lagrange 
equations we find the Langevin equation for (p{x) 

d^df^vix) + d,j d^+'x'x^'ix, x')d'Mx') = J^(x), (16) 

and the Green's function of the system and the medium satisfy the following equations 

^{d^df')G%x - x') = 5ix-x'), 
\d^d'' + Lo^)Glix-x') = S{x-x'). (17) 



In Eq.(16) the noise source and relativistic susceptibility are respectively defined by 



oo 



Jn{x) = -d^ / dujr{uj,x)Y^{x), (18) 







^^^-[x -x')= / dujf^{uj, x)G^{x - x')r{u, x'). (19) 
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Now using Eq.(15) we can find two-point correlation function or Green's function in terms 



of the susceptibility of the medium. We have 

QG 

oo 

X e" 5^ ^ / d"+^x'p{x)G°ix-x')p{x') 

oo 

-:^Id-^+'xJd"+'x' J dcjp^ix)GZix-x')pUx') 

X e . (20) 

Using this recent relation and Eq.Q we find the following expansion for Green's function 

G{x,x') = G\x-x') 

+ [ j (r+h^(r+^z^G\x-z,)[d,,d^,x^^^-{z,,Z2)]G%Z2-x') 



■'p.l^P2' 

In Jn 

7"+l -V . ^_ m+1 m+1 ^ r<0 / 



+ J J J J d^^'z,r+'z2r+'zsr+'z,G'{x-z,)[d^,d^,X^^''%zuZ2)]G'{z2-zs) 
G'izs - z,)[d,,d,,x''''iz3, z,)]G\z, - x') + • • ■ , (21) 



where susceptibility of the medium is given by Eq. ( 19 ) . For a moving homogeneous material, 
that is a medium for which the tensor x^^'^2(x,x') is dependent only on the difference x — x', 
and the coupling function f^{u,x) is independent of x, the Green's function can be written 
in a more compact form in Fourier space as follows 

Gik) = Goik)[l - G,{k)k,,k,,x^^^^{k)]. (22) 
IV. PARTITION FUNCTION 

The partition function of a real scalar field in the presence of a medium is 

S = j D^eT^' = J Dipey'"'"^^, (23) 
where C is given by Eq.Q. The partition function can also be written in frequency [28] 



which is much suitable for our purposes. If we make a Wick's rotation u = iu in frequency 
space, the action will be Euclidean and the free energy is given by = — -InS, where r 
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is the duration of interaction and usually is taken to be quite large. Now from standard 
path-integral techniques we find the free energy at finite temperature T as 

oo 

E = A;bT J]lndet[i^(z/i;x,x')], (25) 

1=0 

where = 2TilkBT/h are Matsubara frequencies and is Boltzman constant. Also 
K^ufj^KjX.') = [z//e(iz/i, x) — V^]5(x — x') and /^(z/;; x, x') = (^^^(iz/;; x, x'). Therefore 

oo 

E = -keT'^trlnlGiiuf,^,^')]. (26) 

1=0 

Now using the expansion Eq.([8]), we obtain the expansion for free energy in terms of the 
susceptibihty function as follows 

E = kBTj2J2-^ / rf"xi---d"x„G°(2z/,;xi-X2)...G°(2z/,;x„-Xi) 

1=0 n=l 

X x{.m,^i)---x{.m,^n), (27) 

where Goiit^, x — x') depends on the dimension of the system. 

V. EXAMPLES 

A. 1+1 Dimension 

Now let us apply the previous formalism to obtain the force induced from a fluctuating 
massless scalar field between two objects with susceptibilities Xi{^) X2{^) in one di- 
mension. In this case the fluctuating field is defined over a (1 + l)-dimensional space-time 
(x = (x, t) G M^^^). In non-relativistic regime, the Green's function of the system and 
reservoir are respectively given by 

-ja;|x-x'| 

G°(a;;x-xO= , (28) 

-ii0'{t-t') 

Gl,{x - x') = 5(x - x')e(t - to , (29) 

and in relativistic regime, these Green's functions are respectively 

G\x,x') = -^\n\x-x'l (30) 
zn 

Gl{x~x')='^Ko{co\x-x'\), (31) 



where Kq{x) = dy 



°° ^r,. cosjxy) 



(x > 0), is the modified Bessel function. Here we restrict 



ourselves to the non-relativistic regime and assume the following form for the dielectric 
function 



eo 



a < X < 6 
6 < X < c 
c < X < (i. 



(32) 



By inserting Eq.(28) into Eq.(27) and considering Eq. (32), up to the first relevant approx- 



imation, the free energy is obtained as 
E = 



^ r r ^g-2i/,|x-x'| 



1=1 



(2^0 



(33) 



now for simplicity we assume that the susceptibilities are position independent, i.e., matter 
has been distributed in Qi and Q2 intervals homogeneously, then 

g-2i/,|x-x'| 



E = -KBTy^x{m)x{m) i dyi j dx.'- 



1=1 



(34) 



for further simplicity we assume that the objects have small dimensions compared to the 
distance between their centers of masses, therefore 

-Ul (d+c—a—b) 



1=1 



(35) 



When the dimensions of objects are not necessarily smaller than the distance between their 
centers of masses we will find the following expression for the free energy 



E 



+ 



-2[d - c) 



where self energies are also included. By defining 



e"^'"''' sinh(2z/ir') sinh(2z/,r")](36) 



a + h = 2ri, h-a = 2r" 

c + d = 2r2, d — c = 2r', r = r2 — ri, 



the force induced by the fiuctuating field can be obtained as 



dE 
Or 



2Xl{yi)X2{yi) -2u,r ■ u/n '\ ■ u/n "\ 

e ' smh(2i/;r ) smh(2z/;r J 



(37) 



(38) 



In zero temperature and in the case that susceptibihties are independent of frequency, the 
force is given by 

F = -^^^(r[0, d-a]+ r[0, c-b]-T[0,d-b]- r[0, c-a]), (39) 
where r[a,2;] is the incomplete gamma function. 



B. 2+1-Dimensional space-time 

In this section we restrict ourselves to a (2 + l)-dimensional space-time (x = (x, t) G M^"*"^). 
In this case and in the non-relativistic regime, the Green's function of the system is given 
by 

Goico, X - x') = '^Mtcu\^ - x'l), (40) 



and the form of the Green's function of the medium is the same as Eq.(29) with the dif- 
ference that now x belongs to a two dimensional space. Now consider two objects with 
susceptibilities xi{^y^) and X2(i^,x), the free energy in first approximation is 

^ = -^E / / d'^'K'oH^-^'\)xiitui,^)x2itu,^'), (41) 

where self energies are ignored. In zero temperature the summation over the positive integer 
/ is replace by an integral according to the rule h |^ — )■ KbTYm^o^ therefore, the free 
energy in zero temperature is 

E = -^, j j j d^^Kl{u\^-^\)xi{tyu^)X2{iy^^). (42) 
In the case that the susceptibilities are independent of the frequency z/, we obtain 

327r'^ J J |x-x'|2 ^ ' 

Using this recent formula we can obtain for example the Casimir energy of two rings |30] . 
Assume that the two ring have radii a and 6, and their centers are separated by a distance 
i?, Fig.(l). This situation can be shown by susceptibilities = Xi^i''^ ~ ^"^^ X2(x) = 

X2<^(^' — b) where r and r' are radial coordinates in cylindrical polar coordinate system. By 



using Eq.(43) the free energy for this geometry is given by 



E^z^i^ r d0 r d0' 1 

327r3 Jq Jq + + b'^ -2aRcose + 2bRcose' -2abcos{e -6')' ^ ' 
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With a simple change in angular coordinates to u = 6 — 6' and v = this expression can 
be integrated to yield the exact closed form [29] 

i?=- ^^^^"^ ^ ^ (45) 

which is the same result obtained from the with weak-coupling regime in scattering formalism 
[30] [3T]. In relativistic regime, the Green's functions of the system and reservoir respectively 
are 



An \x — x'\ ' 



-| p—iui\x—x'\ 

and a similar approach can be followed for obtaining relativistic results. 

C. S+l-Dimensional space-time 

In this section we restrict ourselves to a (3 + l)-dimensional space-time [x = (x, t) G M^"*"^). 
In this case, in non-covariant regime, the Green's function of the system is given by 

1 „-ia;|x-x'| 

G'o(a;,x-x') = -^ -, (47) 

Att |x — x'I 

and the free energy in zero temperature and in the first order of approximation is 

where again we have ignored the self energies which are irrelevant here. If the susceptibilities 
are independent of frequency the expression simplifies further to 

E = -^fd^^[ ^3,. Xi(x)x2(xO _ ( 

647r^ J J |x-x'|3 ^ ^ 

At finite temperature the integral over frequency becomes the Matsubara sum, so the energy 
becomes 

KbT f n /■ ,3 ; , , , coth27rT|x - x'I 
^ = "32^^ J J ^ ^^iW^2(x') ^^-^^^ ^. (50) 

As an example, we consider two spheres, of radius a and b, respectively, with a distance 
between their centers R > a + b. The susceptibilities of the spheres are Xi(x) = Xi5(r — a) 
and X2(x) = X2S{r' — b), where r and r' are radial coordinates in spherical coordinate systems 
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and R lies along the z axis of both coordinate systems, Fig. (2). Then the distance between 
points on the spheres is 

|x - x'l = ^R^ + a2 + 62 _ 2ah cos 7 — 2R{a cos 9 — h cos 9'), 

cos 7 = cos 6* cos 6*' + sin 6* sin 6*' cos (0 — 0'). (51) 



By inserting Eq.(51) into Eq.(49) energy in zero temperature is given by 



E 



XiX2ab 1 - (a - 6)Vi?2 



167ri? l + (a + 6)7i?2' 
which is the same result obtained from the scattering method j31| . 



(52) 



VI. ELECTROMAGNETIC FIELD 



Now we find the Green's function of the electromagnetic field in the presence of some di- 
electrics from which the Casimir force in quite complicated geometries can be obtained 
approximately. For this purpose, the total Lagrangian density can be written in Coulomb 
gauge V ■ A = 0, A° = 0, as follows [32] 

1 1 POO P 

C = -{E' -B') + - duj{Y^{x) - uj^Ylix)) + j duf iu,^)A.Y^. (53) 
Now the interacting generating functional is given by 

POO POO 

+ / dujf{u,i^)AiY^^i + JiAi+ I c/wJ^,iX.,,], (54) 
JO Jo 



where the kernel Kij is defined by 



rj2 

kij = {-^-V'')5,,j-didy (55) 



Now from the well known relation 



and following the same process we did for the scaler case, we obtain the following expansion 
for Green's function in frequency 

G,,(x-x',u;) = G° (x-x',a;) + f ci^G^ (x - zi, a;)[a;2^(a;, zOjC" (z^ - x', u;) + 

Jo. 

/ / d'7.^d^7.2Gl{l^-7.^,u){u^x{^,7.^)]G'l^{7.i-7.2,u)[^^^ , 

Jn Jn 

(57) 
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It can be easily shown that Green's function (57) satisfies the follows equation 

,2 



[-— e(w, x)5ij - V 6ij + didj]Gij{:s., x , u) = 6 (x - x')(5, 



(58) 



where 



+ 1. A similar approach can be followed to find the partition function in 



terms of the susceptibility of the medium as follows 



E = k,Tj:j: 



■1) 



n+l 



n 



rf^xi ■ ■ ■ £y.nG^^^i^{ivu xi - X2)...G'°^i^(iz/i; x„ - Xi 



i=0 n=l 



(59) 



where the free Green's function G°.(x — x', zz/;) satisfies Eq.(58) with e(a;, x) = 1 and uj = ivi. 
By defining r = x — x', we find 



g2 4^j. 



2 



z/;r z/, r 



z/;r vfr^ 3 



(60) 



As an example of the application of Eq.(59), let us find the interaction energy of a system 



composed of two dielectrics with volumes V\ and V2 and the susceptibilities xi and X2, 
respectively. The first relevant nonzero term corresponds to n = 2, therefore 



E = -\kBTS^ f I (i^x(iVG° (x-x',zz/z)G°i(x' -X, zz/z)xi(«z^/,x)x2(«z^«,x') 

Z=0 "^^1 "^^2 



(61) 



Inserting the Green's function (60) into (61), we find 



oo « « 

E = -ksTy^ / rf^X(iVxi(^z^i,x)x2(««^«,x')/i(i/i, |x - x'l), 
1=0 JviJv2 



(62) 



where we have defined 



h{ui, |x - x'l) 



e c 



-|x-x'| ('^l')^ 

■{i 



^/|2 



+ 



2(^)=^ 



6 



+ 



7^}- (63) 



Stt^ '|x — x'l' |x — x'l"" |x — x'p |x — x'l'^ |x — x\ 
In zero temperature, the summation over the positive integer I is replaced by an integral, 
therefore for the case that susceptibilities are independent of frequency, the free energy is 
given by 



E = - 



rf^X / rf3,.Xl(x)X2(xO 



23 



QAn'^ I I |x — x'l''' 



(64) 



For example for two spheres of radii a and b, the distance between their centers R > a + b, 
with susceptibilities Xi(^) = Xi^i^ " and X2(x) = X2'^(^' — b), where r and r' are radial 



coordinates in spherical coordinate systems, the free energy from Eq.(64) is as follows 

23xiX2a^b' 



E 



647r3 



dQ / dQ', 



1 

- V|7' 



X — X' 



(65) 
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Now we use the following formula [?] 

j dndn'\x - x7 = {An fRPPp{d, b), (66) 

where R is the distance between the centers of the two non-overlapping spheres of radii a 
and b, respectively. Here a = a/R and b = b/R, and Pp{a,b) can in general be represented 
by the infinite series 



^^(«' ^) = 1. (2n + 2)! n-p-1) 

= ^ E f ^ M a^("-)6^-. (67) 

2 ^0 \ 2m + 1 y 

We can easily see the following recursion relation holds: 

Pp-ii^^b) = ^^RP+'Ppid,b), P_i = l. (68) 

Using this recent relation, the free energy in zero temperature can be written as 

23xiX2aV 

Eo = ^^P~r, (69) 

where 

-2(a2 - 62)4(^2 + 52 _ 5) ^ (^2 _ 52)2(52(a2 + ^2 _ 44) _ 24(^2 + 52^2) 



10(a4 + (62_i)2_2a2(62 + i)) 



4 



2(5 - 5(a2 + fe2) + 8(a2 + &2)2 _ 4(^2 + 52^3) 
^ 10(a4 + (62_ i)2_2a2(62 + i))4 

In non zero temperature and first order of correction, the free energy is 

23xiX2a'^b^ J. QkBTxiX2a%'^ ,^ ^ , . 

^ ' ^6 (71) 

where 

P_,{a,b) = i[ 1 1 1 + 1 ]. (72) 

4a6 (a-6-l)3 (a + 6 - 1)3 ^a-b+l^ {d + b+l)'^ 

It can be shown that this results are equivalent with the finite temperature weak coupling 
regime in scattering formalism [29l l3T] . 
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A. Covariant formulation 



For a covariant formulation we consider the following Lorentz invariant Lagrangian density 
in Lorentz gauge 



L{x) = -7F^,.F^'" + - / duld^Y^d^'Y^ - cu^Y^] + / dur{co, x)Y^d,A,, (73) 



1 

4-^,-- ■ 2 

where F^''^ = d^A^. — d'^A^ and f^''^{uj^x) is an antisymmetric coupling tensor which couples 
electromagnetic field to the medium. From Euler-Lagrange equations we can show that 
Green's function satisfies the following equation 



UG^,{x -x')- j dS"g^sdXx''''^{x, x")Gp,{x" - x) = g,,5\x - x'), (74) 
where the susceptibility tensor x"^"^{x^x') is defined by 

/■oo 

Xu^.ap{x,x')= / dujr{uj,x)Gl{x-x')r\uj,x'). (75) 
Jo 

For Green's function we find the following expansion in term of the susceptibility tensor 

G,,{x,x') = Gl^{x -x') + J d^z,dh2G%^ix - z^)T''^^^iz,,z,)Gl^iz2 - x') 

J d^z,...d'z,G%^{x - zr)r^^^{zr, Z2)G%^.^{z2 - z,)V'''^\z,, z,)Gl^,{z, - x') + ■ ■ ■ , (76) 

where we have defined V^^^^i^Zi, Z2) = d^^d^^x^^^^'^^^'^i.^i, Z2)- If we are given the susceptibility 
tensor of the medium, we can find the Green's function perturbatively in terms of the 
susceptibility. Having Green's function, we can investigate for example the dynamical energy 
configurations, which is closely related to the dynamical Casimir effects. 



VII. PROCA ELECTROMAGNETIC FIELD 

The Proca equation is a relativistic wave equation for a massive spin-1 particle. The weak 
interaction is transmitted by such kind of vector bosons 133] • The transition from the 
massless electromagnetic field to a field of massive vector bosons is easily done. In contrast 
to the electromagnetic field, the massive spin-1 field will automatically satisfy the Lorentz 
condition dj^A'^ = and the Lorentz condition has become a condition of consistency for the 
Proca field. Therefore the Lagrangian density of the system is given by 

C = --F,,F^^ + ^A,A>^ + dcoiY^ix) - u'yUx)) + j dcof{u;,^)A.Y^. (77) 
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The interacting generating functional can be written in terms of the vector potential and 
the medium fields as 

POO POO 

+ / duf{u,x.)A^Y^^^ + J^A^+ / rfwJa;,^^;^, (78) 
Jo Jo 

where the kernel K^u is defined by 

= [f- V^] V - d,d, - (79) 
Now from the well known relation 

and following the same process we did for the scaler field, we obtain the following expansion 
for Green's function in frequency variable 

^^..(x - x', u) = Gl,{^ - x', uj)+ d^ziGl^ix - zi, u) zi)]G'°,(zi - x', uj) + 

Jn 

/ / rf^zi(i^Z2G'J„(x - zi, cj) [a;^x(^, zi)]G°^(zi - za, u) [u'^xi^, Z2)]G°^(z2 - x', cu) H . 

Jn Jn 

(81) 

Now using Euler-lagrange Equations, we find the equation of motion for the vector potential 
as 

k,,A^ + -j^ —x{co',^)Jdt'e'-'^'-''^—A^{t') = J^ ciu;/(a;,x)yjy,(a;,x), (82) 

where is as a noise or fluctuating field which does not affect the Green's function. The 
Green's function of the above equation satisfies 



[-(— e(a;,x) + -^)V " - d,d,]G^,{x,^ ,00) = 6%x-x')6^,, (83) 



c 



where e(u;,x) = eo[l + x(a;,x)]. It can be easily shown that Green's function (81) satisfies 



in Eq.(83). A similar approach can be followed to find the free energy in terms of the 
susceptibihty of the medium as 

E = ''BTj2J2-^r~ / .2 (^^^ ^1-^2) (^^«;x„-xi) 

1=0 n=l 

X X{m,^i)---X{m,^n), (84) 
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where the free Green's function (^^^^(x — ^'jiui) satisfies in Eq.(83) with e(aj,x) = 1 and 
u = iui. By defining r = x — x', we find 



-Cr 



5^^.(1 + — + 



+ 



+ -6,j'{r), (85) 



where C. = \/ (^)^ — (^)^- an example of apphcations of Eq.(84), let us find the in- 
teraction energy of a system composed of two dielectrics with the volumes V\ and and 
susceptibilities Xi and X2, respectively. The first relevant nonzero term corresponds to n = 2, 
therefore 

E = -\kBT I I c?^x(iVGj;^(x - x , iz/,)G°^(x - x, iui)xi{m, x)x2(«t'/, x'). (86) 

^ "^^1 "^^2 



Inserting the Green's function (p5|) into (86), we find 



E = -keTy^ / c?^xrfVxi(«z^i,x)x2(«t'«,x')/i(z/i, |x-x'|), 

JVi JVa 



(87) 



where we have defined 



g-2C|x-x'| ^4 

h{ui,\x-x'\) = { 
In zero temperature we have 



2C 



5e 



6C 



X — X 



/|2 



X — X' 



/|3 



X — X' 



/|4 



X — X' 



/|5 



+ (88) 



E = -- 



Vl JV2 



rf^xrfVxi(x)x2(x)[(— ) 



+2( 

+5( 



|x — X 

1 



31 
2^13 



31 



h ' |X-X'P^13 



+6( 



+3( 



h |x — X' 



31 

m^i3 



G 



31 



mc 
.mc 



X — X 



X — X 



/|2 1-2 



X — X 



'|2 rs 



-2,0,1 



X — X 



/|2 1-1 



2 2 



X — X 



/|2 



-l,0,i 



31 



n ^|X-X'|6 13 



X — X 



/|2 |0 



.1 1 

2 2 



(89) 



where G^^" (^^l^i' 'b,"^ is the Meijer function [31]. Now for simplicity let us assume that the 
volumes have small dimensions compared to the distance R between their centers of masses, 
such that Iri — ~ -R. Therefore 



E 



CX1X2V1V2 



23 



mc 3 



+ 



,mc 2 3 



+ 



UK' h 16i?6 ' h ' 64i?5 
,mc,n 1 .mc,. 6 



-3 + 27 + 2 log i?) 



+ 



(-1 + 47 + 4 log i?) + 



(90) 



where 7 = 0.577. The first term gives the normal Casimir attractive force and the lowest 
order finite-mass correction tends to reduce the effect. 
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VIII. CONCLUSION 



Quantum field theory in the presence of a medium is investigated for a scalar field and 
massless and massive vector fields in different dimensions. Finite temperature corrections 
are considered and series expansions for corresponding Greens functions and free energies are 
obtained in terms of the susceptibility function of the medium. Some illustrative examples 
arc given showing the applicability and the efficiency of the method. It is shown that the 
Casimir energy in first order approximation in the present method is equivalent with the 
weak coupling regime in the scattering method. The covariant formulation of the problem 
is presented. Finally, the Casimir energy in the fiuctuating Proca field is investigated. 
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Fig.l: Two cylinders of radii a and 6, their centers are separated by a distance R > a + b. 
Fig. 2: Two spheres of radii a and b, their centers are separated by a distance R> a + b. 
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